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Abstract 

The theory of size effects of the properties of nanocrystalline ferroelectric ce- 
ramic or nanoparticle powder allowing for surface tension and depolarization 
field is proposed. Surface tension was included into free energy functional 
and surface energy was expressed via surface tension coefficient. The latter 
was shown to be dependent on temperature due to its relation to dielectric 
permittivity of the nanoparticles. The depolarization field effect was calcu- 
lated in the model taking into account the space charge layer on the surface, 
this space-charge being able to compensate depolarization field in the bulk 
material. 

Euler-Lagrange Equation for inhomogeneous polarization of nanomaterial 
with boundary condition where extrapolation length was shown to be tem- 
perature dependent quantity was solved analytically both in paraelectric and 
ferroelectric phase of size driven phase transition. This phase transition criti- 
cal temperature dependence on the particle size was calculated. Temperature 
and size dependence of nanomaterials polarization and dielectric susceptibil- 
ity was obtained. The possibility to calculate these and other properties by 
minimization of conventional free energy in the form of different power polar- 
ization series, but with the coefficients which depend on particles size, tem- 
perature, contribution of depolarization field and surface tension coefficient 



was demonstrated. These latter effects were shown to influence essentially the 
nanomaterial properties. The comparison with available experimental data is 
performed. 

I. INTRODUCTION 

The characteristic feature of ferroelectric nanomaterials (e.g. the ceramics with nanosize 
grains or powder built by nanosize particles) is known to be inhomogeneous polarization, 
dielectric susceptibility and other properties (see e.g. ref. 0, @, || and ref. therein). The 
properties inhomogeneity is related to the strong influence of the surface, where the proper- 
ties differ from those in the bulk. Because of this surface and correlation energy (polarization 
gradient VP) have to be included when calculating the nanomaterials properties. In phe- 
nomenological approach for ferroelectrics with perofskite structure the additional terms to 
free energy density have the form 5(VP) 2 /2 and <5P S 2 /2A, where P s is polarization on the 
surface, A and 5 are material constants. The variation of free energy functional is usually 
performed to obtain equilibrium polarization. This leads to Euler-Lagrange equation for the 
polarization with boundary conditions, which include the extrapolation length A (see e.g. 
ref. @|) as combined characteristic of surface energy and polarization gradient. Moreover, 
the extrapolation length value reflects the degree of the polarization inhomogeneity. It is 
because the boundary condition has the form of tangent line to the polarization on the sur- 
face (e.g. for spherical particle with radius R the boundary condition is (A^ + P) = 0), 
so that the smaller homogeneity of polarization the larger the extrapolation length and vise 
versa. Therefore it is important to calculate A value rather than to consider it as phe- 
nomenological parameter. While parameter 5 can be estimated via the material correlation 
radius nothing is known about A value. On the other hand calculation of A in Ising model 
includes mainly the change of the average surface coordinate number (see e.g. ref. Q, || 
and ref. therein) that seems to be oversimplification for A value estimation. 

Another problem is related to contribution of depolarization field effect. The influence 



2 



of this field that is able to destroy spontaneous polarization even in bulk materials can 
be decreased by such internal factors as free carriers and domain structure and it can be 
completely compensated under short circuit conditions. It was shown that inhomogeneous 
polarization in the ferroelectric films makes impossible to compensate the depolarization 
field in short circuited films made from single domain ferroelectric, that is perfect insulator 
||. However in the majority of works the influence of depolarization field in nanomaterials 
has not been considered yet. In ref. 0] the influence of depolarization field was considered 
in the model of the multidomain ferroelectric of finite size and space-charge layers on the 
surfaces. However the surface energy of the particles was not taken into account, which 
is known to play the decisive role in the physics of nanomaterials. The important role of 
the surface tension was demonstrated in ||, where observed size dependence of the phonon 
modes frequencies in nanocrystalline PbTi03 was shown to be attributed to the surface 
tension. 

In the present work we propose the model for consideration of depolarization field effect. 
We present the surface energy via surface tension that made it possible to express the 
extrapolation length via surface tension coefficient. This coefficient has been shown to be 
dependent on temperature T, that leads to extrapolation length dependence on T and to 
appearance of some additional temperature dependence of the properties. The ferroelectric 
nanoparticles are supposed to be single domain (e.g. polarization P = P z ) with space- 
charge layer on the surface. The solution of Euler-Lagrange equation for polarization both 
in paraelectric and ferroelectric phases has shown the possibility to obtain the nanomaterials 
properties from free energy expression in the form of power series of polarization like that 
in the bulk, but with renormalized coefficients. These coefficients depend on temperature 
and the particles sizes. Although there is a distribution of these sizes in nanomaterials, so 
that the properties have to be averages with the sizes distribution function, we consider only 
the average radius R of the particles. In many experimental works (see e.g. ref. [[|) 
the properties of nanomaterials were measured as a function of a particle average size. The 
average size can be determined from integrated width of diffraction peak using Scherrer's 
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formula or with the help of TEM method. Both results have been shown to be in a fairly 
good agreement with each other [|TJ. 



II. SURFACE ENERGY AND SURFACE TENSION 
A. Surface energy 

Surface energy is known to be the energy of surface tension on the boundary between 
two media (see e.g. ref. H). It was included into thermodynamic potential written for 



ferroelectric nanomaterial in ref. [|1(| both as a surface energy and as hydrostatic pressure, 
acting on nanoparticles. However these terms related to surface tension were not expressed 
via electric polarization (order parameter of ferroelectrics) , so that it appeared impossible 
to incorporate them directly into Landau free energy that has the form of power series of 
polarization. This incorporation seems to be important to pour light on physical meaning of 
some parameters of thermodynamic theory approach to nanomaterial properties calculation. 

In what follows we will amend this limitation of the previous consideration of surface 
tension contribution to the physical properties of nanomaterials. In general case the surface 
energy can be written in the form |J : 



Esurf = J ^ds (1) 

where /i is the surface tension coefficient and integration is performed over the surface. The 
small change of the surface area ds can be rewritten by the following way 

ds = U xx U yy dxdy (2) 

Here Uu (i = x, y) is a component of strain tensor, the normal to the surface is oriented 
along z axis. Keeping in mind that there is no inversion center near the material surface, 
the strain tensor component can be related to surface polarization P sz via piezoelectric 
coefficients da, namely Uu = du z P z . Substitution of this relation into Eqs. (1), (2) yields: 

E sur f = J [id xxz d yyz P 2 sz dxdy (3) 



Note that Eq. (3) can be used directly for ferroelectric thin films also. Contrary to 
the films in nanomaterials there is another contribution from surface tension related to the 
uniform mechanical pressure p = n/{Ri + R2), where R\ and R2 are the main curvature 
radiuses of the particle surface (i?i,2 — * 00 for the films). For ferroelectrics without inversion 
center this pressure induces polarization via piezoelectric coupling. For ferroelectrics with 
inversion center (e.g. ferroelectrics with perofskite structure) the contribution of the pres- 
sure can be expressed via electrostriction and included into free energy on the base of the 
procedure proposed in []Tl ] . Therefore free energy (with determined stress = —SikP, ex- 
pressed in terms of uniform pressure p = fi/(Ri + R2), polarization P z and strain Uik = 5ikU 
as independent variables) acquires the form: 

5, , 2 3 
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where T c is ferroelectric transition temperature in a bulk material, Qij and C^- are respec- 
tively electrostriction and elastic coefficients in Voigt presentation, E z is the electric field. 
For spherical particles R\ — R2 — R. When taking into account the particle spherical 
symmetry, we assume that the magnitude of oriented along z axis polarization P(r) = P 



depends only upon the radial position r. Analogically to |TTJ and taking into account the 

depolarization field energy, one can obtain the renormalized free energy in the form 

2 



F = 4ir / r 2 dr 



V + b_ pi c_ p sjdpy 

2 4 6 2 I dr J 
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a = a(T — T c ) + 
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Qn + 2Q 



12 



c = 7, 



R C U + 2C 12 ' r ' 6(Cii + 2C 12 )' 

Here E is external electric field, E^ is depolarization field d xxz = d yyz = d is inverse piezo- 
electric effect constant. Hereinafter we will consider ferroelectrics with perofskite structure 
and nanoparticles of spherical form so we will use the Eq. (4b). 
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B. Surface tension 



Let us proceed to the consideration of surface tension coefficient temperature dependence 
related to phase transition from paraelectric to ferroelectric phase. Critical temperature 
dependence of surface tension coefficient when approaching the transition from solid state 
to liquid state is well known M. To the best of our knowledge nobody considered influence 
of structural phase transition like paraelectric-ferroelectric on surface tension. We carried 
out the calculations in the following model. 

In the solids the appearance of polarization manifests itself via the charges on the surface. 
In nanomaterials and in the films the charges can appear in some region near the surface 
where polarization is essentially inhomogeneous. It is known that these charges cause the 
depolarization field in ferroelectrics. In the considered case of single domain nanoparticles 
the only possibility of depolarization field compensation is the space-charge layer on the 
surface that has to include the charges of opposite sign to those related to polarization. The 
possible sources of space-charge layer including free carriers (if any) localized near the surface 
one can find in ref. 0. The influence of the flat surface on the aforementioned charges q can 
be considered by means of the interaction with image charges q* = q(ei — e e )/ (e, + e e ) (see, 
e.g. ref. |I2|), where and e e are respectively the dielectric permittivity of the material of 



particle and surrounding media. All these charges cause the additional energy w s related to 
the surface |§, so that 

A* = »/«.(*,**)«* (5a) 
v 

Here n is the concentration of the charges related to inhomogeneous spontaneous polarization 
(i.e. it is non-zero only in ferroelectric phase), the integration is performed over the volume 
of particle. Taking into account Coulomb interaction between the dipole with moment r d q 
and its image, one obtains the following expression 

nqq f ( 1 



Afi a = -^- / — + )dz (5b) 

P Si J \2z 2z + 2r d 2z + r d ) 



a 
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Here a, r d and / are respectively a lattice constant, minimal distance between positive and 
negative charges and the size of the region of inhomogeneous polarization. Integration in 
Eq. (5) yields 

nq 2 Ei-e e [Aa(a + r d )\ 

This relation can be simplified in supposition r d /a « 1, so that 

nq 2 Ei-E e fr d \ 2 

It should be noted that the consideration we just have carried out can be applied both 
for flat surface and for the surfaces with the principal radii of curvature much larger than 
r d value. This statement follows from the derivation of A/j, s for the case if spherical particle 
performed in Appendixl. 

It should be underlined that since for spherical particles A/x s (T) ~ (£j — e e )/ (ei(ei + E e )) 
the temperature dependence (see below) is the same for all the considered forms. Moreover 
since approximately r d /a < 1, it is possible to say that the expressions (6a) and (Al-6) 
coincide with each other. More generally it would be better to rewrite A/i s in the following 
form: 

Afi s = xnq — — - — r (6c) 

i E e J 

where the value of the coefficient x depends on particle shape, its value being x < 1. 

The surface tension coefficient in ferroelectric phase can represented as the sum of Afi s 
and that of paraelectric phase /z , i-e. 

li = no + A/j, s , (7) 

One can see, that while /i is temperature independent, A/i s depends on temperature 
mainly via £i(T). Really in the considered cases of nanomaterials in the form of nanopowders 
or nanocrystalline ceramics the media surrounding nanoparticles may be considered as non- 
ferroelectric one. In such a case at phase transition temperature €i(T cr ) — > oo and so 



A/i s — > 0, // — ► /io (see (6), (7)). Note, that in the case of two ferroelectric media (e.g. 
BaTi0 3 particles dispersed in SrTi0 3 or PbZr^^Ti^Oa media) both £j and e e depend on T. 
If the transition temperature in surrounding media T = T ec , in these composite materials 
Afi s — > — n q 2 / Ei(rd/ a) 2 at T — > T ec . Note, that the negative sign of Afi s can be valid in 
some temperature region where e e > £j in both aforementioned cases. Keeping in mind 
that surface tension must be positive, it is necessary to add the inequality |A/x s | //xo < 1 to 
Eq.(7). 

Hereinafter we consider only the case of nonferroelectric surrounding media. 

III. DEPOLARIZATION FIELD EFFECT 

We will consider this effect in the model described in previous section. Namely the 
charges localized in the thin surface layer, where polarization is the most inhomogeneous, 
create depolarization field inside the particle. This field can be compensated by space-charge 
layer on the surface (see Section 2). 

The magnitude of depolarization field was calculated on the base of equation written in 
ref. |12] for homogeneously polarized dielectric submerged in a medium with dielectric sus- 
ceptibility e e . To take into account additional charges related to inhomogeneous polarization 
and those from space-charge layer we added the term —Aim a Py a . In the case of external 
field absence this yields 

[(1 - n a )e e + n a \ E da + 4irn a (P a - P Va ) = 0, a = x,y,z (8) 

where Ed is depolarization field, n a is depolarization factor, that depends on the sample 
form and for spherical form n a = 1/3 (see e.g. ref. [0), all the quantities E da , P a , P Va 
being dependent on the coordinates of considered point ~r\ The Py a will be considered as 
variational parameter to minimize the energy related to depolarization field. 
With respect to Eq. (8) this energy can be written as 

« = ^ E / euv = £ E ( (1 _„ a "; + J 2 /(p. - iVJW (9) 
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Equilibrium value of Pvao can be obtained from the condition jjO- = 0, that yields 

Pv*o = ^J dVP^) = P a (10) 

where P is average polarization. 

The Eq. (8) with respect to Eq. (10) gives the following equation for depolarization 
field: 

„ 47m a (P a — P a ) , , 

E da = - .. aK ° , aJ 11 
(1 - n a )e e + n a 

For single domain spherical nanoparticle with polarization along z axis (i.e. a = z, 
n z = 1/3) Eq. (11) gives 

E dz = -a {P z - P z ), a = —— (12) 

1 + 2e e 



Note, that for the film where z coincides with normal to the film surface n z — 1 [Jlz 
so that Eq. (11) leads to the expression obtained in ref. for short circuited film: E dz = 
—Att(P z — P z ). The later is related to the fact that both the short circuit condition and the 
charges on the surface produce similar potential distribution which is able to compensate the 
depolarization field completely (when polarization is homogeneous, i.e. P z = P) or partly 
(when polarization is inhomogeneous) . Note, that Eq. (11) makes it possible to calculate 
the depolarization field for nanoparticles of different form for which depolarization factors 
n a are known. In particular for the particle of cylindrical form with z along the axis of 
rotation n z = 0, i.e. E dz = 0. 



IV. FREE ENERGY FUNCTIONAL AND EULER-LAGRANGE EQUATION FOR 

POLARIZATION 



A. Free energy 



Allowing for the results of previous sections for depolarization field contribution (see Eq. 
(12)), we can write the free energy (4b) in the following form: 
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1 / R 

F = — I 4?r f r 2 ,lr 



o 



-(^P + E )P 
2 4 6 2 \dr V 2 °. 



8vr y fid 2 P 2 rdr 
o > 
(13) 



a = a(T - r c ) + — — - — + a (14) 

Having written Eq. (13) we take into account spherical symmetry and that is why the 
magnitude of oriented along z axis polarization P(r) = P depends only upon the radial po- 
sition r. By virtue of surface tension temperature dependence (see Eqs. (6), (7)) additional 
terms with temperature dependence appear in surface energy (see last term in Eq. (13)) 
and in the coefficient before P 2 (see Eq. (14)). Since \i > and Qn + 2Qi2 > in perof- 
skites both the second term in a (related to surface tension) and the third one (related to 
depolarization field) tends to decrease T c value, i.e. to destroy ferroelectric phase transition. 

B. Euler-Lagrange equation 

Variation of free energy functional (13) yields the following Euler-Lagrange equation for 
static polarization distribution and the boundary condition: 

aP + 6P 3 + cp5 _ s fd*P + 2 dP\ =Eq + aoJ5 
\ dr A r dr J 

[X^ + p) =0 (15b) 

A = ^ < 15C » 

It follows from Eq. (15c), that extrapolation length A depends on surface tension coeffi- 
cient /x(T) given by Eqs. (6), (7), so that A has to be temperature dependent quantity. On 
the other hand fi > in all the materials ||, d 2 > so the sign of A defines by that of 5. 
In what follows we will consider the case 5 > that leads to positive extrapolation length 
A > 0. Note that in the case 5 < one has to add the fourth power of gradient into free 
energy expansion to conserve the system stability. It is seen that it is possible to control 
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extrapolation length magnitude (i.e. boundary condition)by the choice of the materials with 
the smaller or larger piezoelectric constant and surface tension. 



V. SIZE DRIVEN FERROELECTRIC PHASE TRANSITION 



Let us begin with the solution of Eq. (15a) with boundary condition (15b) in paraelectric 
phase where the polarization induced by external electric field can be small enough, so that 
the nonlinear terms in Eq. (15a) can be neglected. In such a case the substitution of 
P — P/r transforms Eqs. (15a), (15b) to the following forms: 



a 6 ~-TT = E o + a oP 

r r dr z 



A 



dP 

dr 



r=R 



r=R 







(16a) 
(16b) 



Keeping in mind, that in the particle center polarization P(0) is finite, one can easily 
obtain the solution of Eq. (16a) and write P(r) as 



P(r) = -(E + a P) 



1 



Rsh(r^a/5) 
7 M(R) 



(17a) 
(17b) 



Averaged polarization can be obtained by the integration of (17a) over particle volume, 



R 



i.e. P = 3/R 3 j P(r)r 2 dr, that yields 



P(R) 
W(R) 



E W(R) 
a - a W{R) 

i 3 5 R^a7Sch(R^aj6) - sh(R^/S) 



R 2 a M(R) 
Substituting Eq. (18a) into (17a), one obtains 

E 



(18a) 
(18b) 



P r) 



, R sh(r/a/5) 
r M(R) 



(19a) 



a - a W(R) 

It is clearly seen that polarization in paraelectric phase is proportional to external field, 
the coefficient of proportionality being the inhomogeneous dielectric susceptibility 
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XPE[r) 



a - a W(R) 



Rsh(r^a/5) 
7 M(R) 



(19b) 



It is follows from Eqs. (18a,b) and (14) that for bulk material (R — > oo) W(R) — > 1 and 
P{R) -> E /(a{T - T c )) as it has to be. 

The average value of susceptibility can be extracted from Eq. (18a) as the coefficient 
before Eq. It follows from Eq. (18a) that the susceptibility becomes infinitely large at the 
condition: 

a-a W(R) = (19c) 

Allowing for that both a and W(R) depend on particle radius R and temperature T, 
the condition (19c) can be satisfied at some critical radius R cr (temperature is fixed) or 
at some critical temperature T cr (radius is fixed). These critical values correspond to size 
driven ferroelectric phase transition in nanomaterials. Its consideration without contribution 



of depolarization field and surface tension was carried out earlier in e.g. ref. [§[], 0, |L3 
|14| . Comparative analysis of the results with and without these factors contribution will 
be performed later. 

Because of complex enough form of W(R), that depends also on T via temperature 
dependence of parameter a (see Eq. (14)) it appeared cumbersome to calculate T cr and 
R cr analytically in general case. Because of this we performed the calculation numerically 
introducing the following dimensionless parameters and functions: 

^ = VZ' 0O=-rr> Po = —ET \ n , ' ( 2 °) 



4vr' u 4tt ' ru ttR \C 11 +2C 12 

VAtt6 xnq 2 e< - s e 

A = — rr, V = , An FE (T,R) 



Pod 2 ' /'o ' " Si{£i + e e )' 

The results of T cr calculations is represented in Fig. 1. 

The increase of fiQ value at other fixed parameters leads to the decrease Aq, that results 
in critical temperature decrease and in critical radius increase (compare the curves 1 and 
2 and their crossing points with abscissa axis). Also one can see, that the decrease of 
surrounding medium dielectric permittivity e e leads to R cr decrease and T cr increase, these 
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effects being larger for smaller A value. Dashed-dotted curves correspond to approximate 
analytical formulas obtained by the following way. 

Estimations have shown, that for e e < 100 and R > 10 nm the contribution of ao into 



a value is the largest one, i.e. further we can put a ~ ao- Keeping in mind that \j5/a can 
be estimated as correlation radius value, that has to decrease due to depolarization field 



contribution, it is reasonable to suppose that y5/a value is about several lattice constants. 
Therefore for R^Ja /5 > 3 5 one can derive from (19c), allowing for that at transition 
point /i — > /io, the following approximate analytical expression for the critical temperature 
of size driven phase transition: 

T cr (R) « T c (l - § - (21a) 



R R — R c 



Here we introduce the parameters 



i?oA V47T Ro D Ro f 3 v /a ~ 

tic — - — — 7=, Ml — Po~pr, tC\ — 



i? c can be considered as correlation radius renormalized by contribution of depolarization 
field and surface tension. i? L is correlation radius renormalized only by surface tension. 

For arbitrary T < T c approximate formula for critical radius R cr (T) can be obtained 
from Eq. (21a), as the solution of the quadratic equation with respect to T cr = T,_R = R cr . 
The analytical expression can be obtained for the some cases. 



In particular, for the large enough particles when R^a /5 » 1 and Rj R c » 1, one can 
neglect R c in the denominator of the third term of Eq. (21a) and formula for T cr (R), _R cr (T) 
dependences can be written as 

T cr (R) « T c (l - =T e (l- ^jp- (1 - T/T c )) , (21b) 



Rcr{T) « — Po + 



6 (1-T/T C ) V Aov^+v 7 ^; ' 
If R L « R x and i? A / (1 — T/T c ) « R c approximate formula for T cr (R), R cr (T) depen- 
dences can be written as: 

T cr (R)~T c (l- \ ), R cr ^R c (21c) 

V it — tier' 
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Thus, in this particular temperature range critical radius does not depend on temperature. 

It is seen from (21b) that (T cr — T c ) is inversely proportional to the particle size and 
e.g. for ferroelectrics with perofskite structure T cr < T c because Qu + 2Qi 2 > (although 
Qu < 0). It should be noted that the value of (T cr — T c ) is proportional to the coefficient 
of surface tension, the first and the second term in the brackets being related to surface 
tension and surface energy respectively. One can see, that Eq. (21b) describes the main 
features of T cr behavior: T cr = T c at yU = 0, T cr decreases with /i increases and T cr value 
increases with the particle size increases (T cr — > T c at R — > oo). The results of more detailed 
calculation of T cr (R) dependence given by Eq. (21b) is shown in Fig.l by dashed-dotted 
line. One can see that the coincidence between calculations on the base of Eq. (19c) and 
on the base of approximate formula (21b) is very good for A = 250, for A = 25 the T cr (R) 
dependences given by curve 1 and dash-dotted line are similar, however quantitative values 
of T cr for both curves are different. 

The equation (21b) at T = does correspond to the crossing points of the curves in Fig. 
1 with abscissa axis. But it is possible to obtain R cr at any fixed arbitrary temperature with 
the help of Eq. (21b). 

It follows from Eqs. (21) that critical parameters of size driven ferroelectric phase tran- 
sition in nanomaterials are defined completely by physical characteristics of the material 
such as Curie- Weiss constant Cw = 4:7i/a, electrostriction constants Qij, elastic modulus 
Cij, piezoelectric constant d lz = d and surface tension coefficient /i . This appeared possi- 
ble due to the fact that extrapolation length that usually is a phenomenological parameter 
was expressed via /i and d coefficients and surface tension was related to polarization via 
electrostriction effect (see section 2). 
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VI. PROPERTIES IN FERROELECTRIC PHASE AND FREE ENERGY WITH 

RENORMALIZED COEFFICIENTS 



In ferroelectric phase R > R cr , T < T cr the nonlinear term in Eqs. (13), (15) can not be 
neglected. The simplest way to take it into account is to look for polarization on the base 
of direct variational method. We choose solution (19a) as a trial function and amplitude 
factor will be treated as variational parameter, therefore 



P FE = p v | 1 _ 



Rsh(rJa/5) 



(22) 



r M(R) 

Here Py is variational parameter that represents the amplitude of polarization space distri- 
bution. 

After substitution of the trial function (22) into integral (13) one can obtain the free 
energy in the following form: 



P - Ar P 2 _l Br p4 , Cr p6 p p 



(23a) 



Here 



R 



r 2 dr 



a 1 



Rsh(ryja/S) 
7 M{R) 



d_ R shjr^Ja/S ) 
+ 1 dr 1 /• M(R) 



-a Q (W(R)) 2 + 



5 3 / sh(R^a/5) 



XR 



M(R) 



R 



B^ h llr^r[l- RSh{rH5) 



R 3 



E R = E W(R), 



r M(R) 



Cr ~ R 3 ( dr y r M(R) 



(23b) 



where parameters a, b and c are given by Eq. (14). 

Keeping in mind, that average polarization is P = P V W(R) ) we can rewrite Eq. (23a) 
as follows 



A R P 2 B r P 4 C r P 6 
p _ n I tl I ii pp 

W 2 {R) 2 W\R) 4 W e (R) 6 
15 



(23c) 



It is worth to underline that the terms in Eq. (13) which correspond to polarization gra- 
dient, depolarization field and surface energy contribute to A R (see respectively the second, 
the third and the forth terms in Ar). Therefore the main peculiarities of nanomaterials 
properties is expected to be related to the first term in Eqs. (23a), (23c). 

It is obvious that Eqs. (23a) and (23c) have the from of power series of the amplitude 
of the polarization space distribution and average polarization respectively. In general case 
all the coefficients (23b) depend on a particle size and temperature. 

The main advantage of the Eqs. (23a, c), that they make it possible to obtain P v and 
P by conventional minimization procedure. In particular for the phase transitions of the 
second order 



and the polarization profile can be easily obtained with the help of Eq. (22). 

The physical properties e.g. dielectric susceptibility x( r )> X an d pyrocoefficient II(r), II 
can be calculated as the derivative of polarization over E and T respectively. 

The results of calculations of average dielectric susceptibility dependence on the particle 
size and temperature are depicted in Figs. 2 and 3 respectively both for paraelectric (R < 
R cr , T > T cr ) and ferroelectric (R > R cr , T < T cr ) phases for different values of dimensionless 
extrapolation length A . It is seen that inverse susceptibility linearly depends on T like 
that in Curie- Weiss (C-W) law (see Fig. 3), although in ferroelectric phase there is small 
deviation from linearity for small extrapolation length (see curves 1). In paraelectric phase 
1/x ~ (R cr / R — 1) while in ferroelectric phase the dependence of x on R 1 declines from 
linearity more strongly, especially for /i(T)//i — 1 = r]Afi FE (T) ^ (see solid curve 2), 
although general view resembles that of C-W law (see inset to Fig. 2) with maximum at 
R = R cr . It follows from Figs. 2,3 that at /i = x = Xb- (because at fi — > A — > oo, that 
corresponds to the bulk) and the decrease of i] value makes the dependences linear, i.e. 1/x ~ 
(R cr /R — 1) and 1/x ~ (T — T c ) (see dashed curves in the Figs.). Thus all deviation from 
linearity in ferroelectric phase is related to A/i FE (T) contribution. The influence of dielectric 




(24) 
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permittivity e e of the surrounding medium appeared to be strong enough (compare solid 
and dotted curves in Figs. 2,3). This effect is related to the contribution of depolarization 
field and temperature dependent part of surface tension. 

Practically the same type of influence of /xo, A/i and e e on the averaged spontaneous 
polarization one can see from Figs. 4 and 5, where the temperature and size dependence of 
averaged polarization are depicted, the crossing points of the curves with abscissa axis in 
Figs. 4 and 5 being respectively R /R cr and T cr /T c values. One can see, that the dependence 
of this quantities on extrapolation length is a good agreement with approximate formulas 
(21b, c), where the second term in the brackets reflects the influence of extrapolation length 
i.e. the larger A the larger T cr and the smaller R cr value has to be. This behavior is related 
to the fact, that larger extrapolation length corresponds to more homogeneous polarization 
like that in bulk material, so that the closeness of T cr to T c has to be expected. It is evident, 
that the destruction of such "strong" polarization can be achieved for small enough particles, 
i.e. for small R cr value. It is worth to underline, that the calculation of T cr and R cr with 
the help of Eqs. (21b, c) leads to the following values: T cr = 0.44T C , R cr = 0.036i? (for 
A = 25), T cr = 0.755T C , R cr = 0.0785^0 (for A = 250), which are cloth to those obtained 
from crossing points in Figs. 4, 5. This speaks in favor of statement that the accuracy of 
formulas (21b, c) is not bad at all. 

The calculations of space distribution of polarization (see Eqs. (22), (24)) had shown, 
that the contribution of temperature dependent part of surface tension can be large enough 
(compare solid and dashed lines in Fig. 6). It is seen, that r]Afj, FE (T) ^ flattens the curves 
and decreases polarization, these two effects can be related to increase of extrapolation length 
and surface tension respectively with fJ>(T) //io = 1 +r]AfiFE increase. Therefore the influence 
of surface tension and its change under externally controllable parameters (e.g. temperature) 
can manifest itself mainly in space distribution of physical quantities, the smaller the particle 
or extrapolation length the larger the effect can be (see Fig. 6). 

The considered temperature and size dependences of properties were calculated numer- 
ically because of the complex form of the coefficients (23b). These dependences can be 
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more clearly seen after analytical calculation of integrals (23b). In particular for A R the 
integration yields : 

A R =[a-a W(R)]W(R). (25) 

It is seen that expression in the brackets completely coincides with left hand side of Eq. 
(19c), so that at T = T cr or R = R cr A R = 0, Py = 0, P = as it has to be. 

The expressions for B R and C R appeared to be very complex (see Appendix 2). But 
all the coefficients can be simplified allowing for the estimations made in previous section 
for T cr and R cr approximate calculations. Since the accuracy of these formulas was shown 
to be good, it is worth to write the approximate formulas for the coefficients A R , B R and 
then for the physical properties. Allowing for the approximate expression for W(R) at 



RJclq/S > 3 -J- 5 one can rewrite A R approximately in the following form: 



A R ^a(T-T c ) + ^ 



4 gn + 2Q 12 3 d?y/5a^ 



Cn + 2Ci2 d 2 fj, + a/5oo 
Substitution of Eqs. (21b) for T cr or R cr into Eq. (26a) yields 



(26a) 



A R = a (t - T cr + ^(T c - T cr )^) , (26b) 
A R = a(T-T c )(l-^±^^). (26c) 

One can see, that for the case when A(j, s //i <^ 1 expression (26b) corresponds to that 
for C-W law, and (26c) contains the simple multiplier (1 — R cr /R). But in general case 
situation is more complicated. Indeed, because A/i s depends on the material dielectric 
permittivity (see Eq. (6a)) it can be expressed via A R . In particular it is possible when 
(£j — e e )/(£j(£j +£ e )) ~ that is proportional to A R . In such a case (26b) and (26c) can 
be represented as 

A R = a{T - T cr )A R1 (27a) 



A R = a(T-T c )[l-^)A R2 (27b) 
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Ari l-m7(T c -T cr )' Am l + ap-T^Rcr/R (2?C) 
Allowing for that e.g. for BaTi0 3 rj 10 2 and a ~ 10~ 5 one can see, that A R1 1 and 

A R2 ~ 1 although the accuracy of this estimation for A R1 is better than for A R2 . In what 

follows we will use Eqs. (27a), (27b) with A Ri = A R2 = 1- 

The formulas (27a) and (27b) can be applied respectively to study the temperature 
dependence for some fixed radius of the particles and to study the size dependence at some 
arbitrary fixed temperature. 

Keeping in mind, that in adopted approximation a ~ ao and so it is independent on 
temperature, the coefficient B R also independent on T (see Eq. (23b) and Appendix 2), so 
that the properties temperature dependence defines completely by (27a). Moreover, because 
W(R) m 1, so P v m P (see Eq. (24)) and E R « E (see Eq. (23b)). In such a case the 
susceptibility has the conventional form 

X(T)= 1 T>T cr (28a) 

a(T - T cr ) 

and 

X(T)= 1 T<T cr (28b) 

2a(i cr - i J 

respectively for paraelectric and ferroelectric phase of size driven transition. The tempera- 
ture dependences given by Eqs. (27a), (28a), (28b) correspond completely to those in bulk 
materials but with substitution T cr for T c . These dependences are shown in the Fig. 3 
by dash-dotted curves. One can see, that the approximate formulas describe the accurate 
curves good enough, especially curves numbered by 2, where dash-dotted curve completely 
coincides with solid curve, although the coincidence is not so good for the curves 1. Tak- 
ing into account, that for ferroelectrics E{ ~ 47rx, it is easy to obtain after substitution of 
T cr (R,R cr ) from Eqs._(21b) into Eq._(28), that at fixed temperature T the dependence of 
nanoparticle dielectric permittivity over particle radius R can be represented as: 

*< r -* ) = 2(l-£(7W H>H " in <28C) 
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Here e b = C W /(T C - T). 



The spontaneous polarization approximate formulas can be obtained from Eq. (24) at 
W(R) ~ 1 and Br (3 (see Appendix 2). This yields for average polarization 



This formula completely coincides with that for bulk material when substituting T cr for 
T c . The approximate expression (29) describes exact curves P(T) not bad at all (compare 
solid and dash-dotted curves in Fig. 5 which completely coincide for A = 250). Note, that 
these curves describe also the temperature dependence of inverse pyrocoemcient II because 

— dP a 1 , . 

n = 5r = -^ < 30 > 

For the phase transitions of the first order the spontaneous polarization, the dielectric 
susceptibility and the pyrocoefficient depend also on Cr coefficient given by (23b). It was 
shown (see Appendix 2) that Cr ~ c, so that all the physical quantities for nanomaterials 
can be described by the same formulas as for bulk materials, but with the coefficient before 
P 2 in the form of (27a) or (27b). The substitution of Eq. (27b) will make it possible to 
describe approximately size dependence of the properties. The description on the base of 
the approximate formulas appeared to be good enough especially for the curves 2 (compare 
solid and dash-dotted curves in Figs. 2, 4). 

The space distribution of the properties in nanomaterial can be calculated with the help 
of Eq. (22). Since in adopted approximation a m a is temperature independent quantity 
and P = Py the properties profiles will be described by the expression in the brackets of 
Eq. (22), i.e. their form will be the same as that depicted in Fig. 6. 

VII. THE DEPOLARIZATION FIELD AND SURFACE TENSION COEFFICIENT 
SIZE AND TEMPERATURE DEPENDENCIES 




(29) 



In this section we will consider in details size and temperature effects in depolarization 
field and surface tension coefficient. 
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A. The depolarization field 



It was shown in Sect. 3, that for nanoparticles with spherical form the depolarization 
field is given by Eq. (12) in the model, when it can be partly compensated by space-charge 
layer. If dielectric permittivity of surrounding media e e independent on T, temperature 
dependence of E dz = E d will be related to that of (— P z + P z ) = (— P + P). Allowing 
for space distribution P(r) (see Eq. (22)) one can calculate E d (r) distribution, which is 
depicted in the inset to Fig. 6. It is seen that for infinitely large extrapolation length 
(A ~ l//iO) A*o — 0) that corresponds to practically homogeneous polarization, E d = 0. For 
smaller extrapolation length values when polarization is inhomogeneous E d < inside the 
particle and transforms into E d > near the surface. The absolute values of E d decrease 
with A increase (compare the curves 1 and 2 in the inset to Fig. 6). Therefore depolarization 
field tends to decrease the polarization inside the particle and to increase it near the surface. 
Since in the considered case of positive extrapolation length polarization near the surface 
is smaller than inside the particle, depolarization field tends to flatten the polarization in 
nanoparticle, i.e. to make it more homogeneous. 



B. Surface tension coefficient 



Let us discuss now the size and temperature dependence of surface tension coefficient in 
ferroelectric nanomaterials. In the adopted model surface tension coefficient is given by Eqs. 
(7), (6), so that its size and temperature dependences are defined by that of nanoparticle 
dielectric permittivity Si in supposition that permittivity of the surrounding media e e is 
temperature independent quantity. Because of Si dependence on the temperature T and 
particle radius R the quantities A/i s and Afi FE from Eqs. (6), (20c) have to depend on 
these parameters. One can see the dependence A[i FE (T) in Fig. 7, and A/j, FE (R) in Fig. 8 
for the several values of e e . Negative sign of Afi FE corresponds to the case £j < e e) that 
realizes at e e > 300 both for Aq = 25 and 250 (see inset). The existence of the maximum in 
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curve 3 is related to the temperature region where £j and e e is close to one another. Positive 
sign of AfipE corresponds to the case > e e (see curves 1, 2), the curve 1 being very close 
to the respective curve for 1/x dependences on temperature (see Fig.3) and on the particle 
radius (see inset to Fig. 8). 

One can see from Figs. 7, 8 that \A/j,fe\ ~ 10~ 3 -f- 10~ 2 , i.e. the relative contribution of 
A// s to the surface tension \x is equals A// s //i = ^A/ip^lO -1 4- 1, for 77 = 100 that we used 
when plotting all the Figures 1-6. 

Therefore for nanomaterials A/j, s (T) can be both positive and negative quantity, its abso- 
lute value \A/j, s (T)\ < fiQ. This different possibilities depend on the material characteristics 
(particles size, extrapolation length, transition temperature in bulk) and dielectric constant 
of the surrounding media. 

The temperature dependence of surface tension coefficient results in extrapolation radius 
temperature dependence because A ~ 1/ fi (see Eq. (15c)). Since in paraelectric phase \i = //q 
is temperature independent quantity, we depicted X/Xpe = 1/(1 +Afi a /jj, ) in Fig. 9 as the 
function of temperature. Note, that Xpe = RqXq- One can see, that in paraelectric phase 
(T > T cr ) X/Xpe = 1, while in ferroelectric phase (T < T cr ) its behavior depends strongly 
on e e value as it has been already discussed above, allowing for A// s /// = r]Afi FE . The 
same type behavior was obtained for X/Xpe dependence on particle size with X/Xpe = 1 at 
R < R cr (see inset 1 to Fig. 9). Note, that in the case A/x s > the extrapolation length 
value decreases. Under the condition when A/x s < extrapolation length has to increase 
resulting into more homogeneous polarization in the nanoparticle. 

VIII. DISCUSSION 

1. Let us compare our theoretical calculations with experimental dependences T cr (R) 
and el(R) for BaTi03 and PbTiOs nanomaterials [0, ||16|| , 0. 

Taking into account that R\ ~ l/\/l + 2e e (see Eqs. (12), (20), (21a)), i.e. decreases 
when e e increases, and R c saturates when e e increases, one can obtain that two cases are 
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possible. For ceramics with e e » 1, one can suppose that R\ tC << Rl and use Eqs. (21b), 
(28c). In contrast, for powder samples e e ~ 1, one can suppose that R\ » Rl and therefore 
Eqs. (21c) is valid. Really, the estimation of the parameters for aforementioned materials 
(see, e.g. JHJ, JTT], |Tj|]) have shown that Eqs. (21b) and (21c) can be used for the ceramic 
and powders description respectively. 

Below we regard critical radius as the fitting parameter, so Eqs. (28c) depend on the one 
fitting parameter R cr = R cr (T) at fixed room temperature T. The comparison of calculated 
on the base of Eqs. (28c) with the experimental £l(R) dependence for BaTiOs ceramics 
is represented in Fig. 10. It is seen good agreement between experimental and calculated 
data. The difference between obtained critical radii value can be the result of difference 
in intergrain media due to different technological processes of BaTiOs ceramics production 
(see 0, and ref. therein) which can result in difference of surface tension coefficients. 

Critical temperature T cr (R) depends on the two parameters: radii R cr and R\. The com- 
parison of (21c) with the experimental dependence T cr (R) for powder BaTiOs and PbTiC>3 
samples is represented in Fig. 11. One can see that the theory fits the experimental data 
rather well. 

Using the fitting parameters values, given in the captions to figures 10 and 11, the 
known values of the electrostriction, elastic stifness, Curie-Weiss constants and criti- 



cal temperatures (see, e.g. Hl7[ , [fUJl ) and the following values of the polarization gra- 
dient coefficients for BaTi0 3 and PbTi0 3 in CGSE-units (see, e.g. fL8"f ): <5(BaTi0 3 ) = 
5 • 10~ 15 cm 2 , <5(PbTi03) = 5 • 10~ 16 cm 2 ,one can obtain from Eqs. (15c), (21a), the the 
following data for physical quantities used in the theory: /i = 3.11 • 10 5 din/cm 2 , Xpe = 
16 nm (BaTiOs ceramics with e e = 100 and d = 10~ 7 CGSE-units); /io = 3.28 • 
10 5 din/cm 2 , Xpe = 50 nm (BaTiOs powder with e e = 1 and d = 0.55 ■ 10~ 7 CGSE-units). 
Using the known value of the surface tension of PbTi03 powder fi = 5 ■ 10 4 din/ cm 2 || 
and e e = 1 one can obtain d = 1.20 ■ 10~ 7 CGSE units, Xpe = 6.9 nm. 

One can see that the obtained quantities have reasonable values known for these param- 
eters (see, e.g. fl5H , Hl6| , fL8|). 

23 



2. The surface tension and depolarization field influence were taken into account in 
the thermodynamic approach to the investigation of ferroelectric nanomaterials properties. 
Surface energy was expressed via surface tension coefficient that results into extrapolation 
length dependence on this coefficient. Its value was obtained in the framework of the model 
taking into account the space-charge layer on the particle surface and shown to be tem- 
perature dependent quantity. In the same model the depolarization field was taking into 
consideration. This field was shown to be partially compensated by the surface charges and 
proportional to the deviation of the polarization from its mean value. Coefficient of pro- 
portionality is defined by the nanoparticle shape and surroundings dielectric permittivity. 
This lead to the essential difference between dielectric properties of the ceramic and powder 
samples because of the surrounding media properties difference. 

The study of the physical properties of the ferroelectric nanomaterials such as transition 
temperature, critical particle size, spontaneous polarization and dielectric susceptibility were 
performed on the base of Euler-Lagrange equation for the polarization with the boundary 
condition with the temperature dependent extrapolation length of the positive sign. It 
should be underlined that its negative sign can arise from the negative sign of the coeffient 
before square polrization gradient. But in this case the higher power gradient terms has to 
be included into the free energy in order to stabilize the system. 

The possibility to calculate physical properties of the nanomaterials by the minimiza- 
tion of the "bulk" free energy density with renormilized coefficients, which depend on the 
particle size, temperature, surface tension and depolarization field characteristics has been 
demonstrated. The developed theory gives the basis to the empirical expression proposed in 



15] for the BaTi03 and PbTi03 transition temperature dependences on the particles size 



and fits rather well the available experimental data. 
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IX. APPENDIX 1 



Let us consider the point charge q located at the distance r from the center of the 
dielectric sphere with radius R > r and dielectric permittivity The outer space is also 
regarded as dielectric with permittivity e e . 

In order to find the electric potential (p l q ' e of this system one have to solve Laplace 
equation with the appropriate boundary conditions. In this case it has the following form: 

A V?r = °> e&tp* =eeVtf D (Al-1) 



r=R 



r=R 



r=R 



Having solved this equation in the spherical coordinates it is easy to derive the surface 
contribution to the potential ip q : 

<(ro , r , 0) = £ ~ ~ ~ + - ^ (cos W ) ; r<* (Al-2) 

where (r, 0, </?) are spherical coordinates and P m (cos(#)) is Legendre polynomial of m-th 
order. This potential is determined by the surface bound charges induced by charge q. Thus 
the additional surface energy can be calculated as interaction energy between charge q and 
induced charges. Considering the electric dipole with the dipole moment q r d directed along 
polar axes, one can write the interaction energy with respect to P m (cos (9 = 0)) = 1 as 
follows 

w s (r ) = q (i<P s q (r , r , 0) - ip s q (r + r d , r , 0)) - (<p*(r , r + r d , 0) - tp s q (r + r d , r + r d , 0))) « 

^ / ^2 £ i ~ £ e m + 1 m 2 7"o m , . 

~ [qTd) ~h^Bi + {m + l)E e (A1 " 3) 

Here we used the strong inequality r d « r . Then the surface energy can be represented 

in the form of the equations (5a) (5b) 

R-a 



Having performed the integration, one can obtain 

A// = n ^£^££ £ - + 1 - ^ (l - ^-) 2m+1 (Al-4) 

R 2 e t ^ me i + (m+l)s e 2m+l\ Rj K ' 
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The majorant series of the series from (Al-4) is equal to 

^tfr + (Al-5) 
Here b = (1 — a/R) 2 . Sum of the series (Al-5) can be find in the following way 



Y (m + i)b m =Y ^1 = — (y f+A =±(J-) = — L_ 

h { } h db db\h ) db{i-b) {i-by 

Using this relation it is easy to find the evident dependence of the series (Al-5) on a/R ratio. 
Owing to the inequality a « R we finally get the following expression for the additional 
surface tension (Al-4) 

Afl =n 8\7 g ) £i{£i + £ e y (A1 " 6) 



X. APPENDIX 2 



Let us proceed to the integration in the equations (23b). Owing to the boundary condition 
(15b) integration by parts of the second term in A R leads to the following equation for the 
coefficient A R 

« R 



sinh(ry a/ 5) 



+5 



A R + a om R)f=-fr* d r(l ^ ^ 

3f ' 



a 1 



sinh(ry a/5) 
M(R) r/R 



+ 



R d 2 sinh(rya/<5) 
~dr^ M(R) 



r dr 1 — 



sinh (rya/5) 
M(R) r/R 



(A2-1) 



Since this integral is equal to W(R) one can get the expression (25) for the coefficient A R . 



For the sufficiently large particles radius R » yd /a it is easy to substitute cosh(:r) and 
sinh(x) for exp(x)/2 and to rewrite expression for W(R) as follows: 

3 



W(R) 



1 - 



Z(t + i)-t 



(A2-2) 



Here 



f = Rja/S, i = A 1 
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Integration in the expression for Br from the equation (23b) leads to the following 
coefficient form 

B " 1 " o„»Lu + m 2 M{Rf + 96^ 3 (Ccosh(0 - sinh(O) - 



b 8eM(R) 4 

-12£ 2 cosh(2f ) + 3£ 2 cosh(4£) - 36£ M(R) 2 sinh(2f ) - 

-72£ 2 M(R) Sih(0 + 24£ 3 Sih(2£) + 24£ 2 M(i?) Sih(3f ) - 12£ 3 Sih(4£)) (A2-3) 

Here and Sih(£) = Jj; d( smh(() / ( is the integral hyperbolic sine function. When using the 
strong inequality (21c) in the form £ >> 1 and the asymptotic form of Sih(£) ~ exp(£)/(2£) 
one can suppose that M(R) = exp (£)(£. + l)/2 and obtain the following expression for B R 

B «™ h { l -f(^ {U+l) ) (A2_4) 
It is seen that under the condition 

\*[% = l»\ (A2-5) 
V o 

coefficient B R is close to b. 

Coefficient Cr also can be expressed in terms of the hyperbolic functions and the integral 
hyperbolic sine function, but the final expression is very cumbersome and omitted for the 
sake of briefness. Using inequalities (21c) and (A2-5) we can suppose that 



Rsinh(rya/5) R2sinh(rya/5) 
1 ~ 7 M(R) ~ 1 ~ 7 exp(£) {£ + 1) ~ 1 

Therefore when using the explicit form (23b) of the coefficients Br and Cr their approximate 
expression can be easily simplified to 

B R « b 

Cr^C 

The first of these approximate equalities coincides with (A2-4) owing to strong inequality 
(A2-4). 
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FIGURE CAPTIONS 



Figure 1. The size dependence of transition temperature at 6 = 0.005, e e = 95 and 
p = 0.05 for different A and rj: A = 25 (curves 1),A = 250 (curves 2), rj — 100 (solid 
curves), i] = (dashed curves), e e = (dotted curves), approximate formulas (dash-dotted 
curves), po = (dash-double-dotted curves). Inset: the same dependences over the inverse 
radius. 

Figure 2. The dependences of the inverse static susceptibility on inverse radius (xb is 
bulk static susceptibility) at fixed temperature 6 (T/T C — 1) = —0.005. Inset: the same 
dependences of x/Xb on R/Rq. The parameters for the different curves correspond to those 
in Fig.l. 

Figure 3. The temperature dependences of inverse and direct (inset) static susceptibility 
at fixed radius R/Rq = 50. The parameters for the different curves correspond to those in 
Fig.l. 

Figure 4. The size dependence of spontaneous polarization on the inverse and direct 
(inset) particle radius (Psb is bulk polarization) at 6 (T/T C — 1) = —0.005. The parameters 
for the different curves correspond to those in Fig.l. 

Figure 5. The temperature dependences of spontaneous polarization over T/T c (Pg b is 
bulk polarization at zero temperature) at fixed radius R/Rq = 50 and 9 (T/T C — 1) = —0.005. 
The parameters for the different curves correspond to those in Fig.l. 

Figure 6. The spatial distribution of the spontaneous polarization and the depolariza- 
tion field (inset) at fixed temperature 9 (T/T C — 1) = —0.005, R/Rq = 40. The parameters 
for the different curves correspond to those in Fig.l. 

Figure 7. The temperature dependence of surface tension coefficient at fixed radius 
R/Rq = 50, 9q = 0.005, A = 25 and p = 0.05, rj = 100 and different e e values: e e —l 
(curves 1), e e =100 (curves 2), e e =300 (curves 3), e e — > oo (curves 4). Inset: the same 
dependence for A = 250. 

Figure 8. The size dependence of surface tension coefficient Afx FE , parameters are the 
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same as in basic plot in Fig. 7. Inset: the same dependences for 1/x in FE-phase. 

Figure 9. The temperature dependence of extrapolation length (basic plot), parameters 
are the same as at basic plot in Fig. 7. Inset: the size dependence of extrapolation length at 
fixed temperature 6q(T/Tc — 1) = —0.005. 

Figure 10. The comparison of experimental (open squares 0, triangles ||) and theoret- 
ical given by Eq. (28c) (solid curves) dependences of the static permittivity (basic plot) and 
impedance (inset) over particle diameter D for BaTiC>3 ceramics. The fitting parameters: 
Eb = 1450, D cr = 700 nm (basic plot), D cr = 1500 nm (inset). 



Figure 11. The comparison of experimental (open squares ||15|| , fllql ) and theoretical 
given by Eq. (21c) (solid curves) dependences of transition temperature Tqr over parti- 
cle diameter D for powder BaTiC>3 (basic plot) and PbTiC>3 (inset) samples. The fitting 
parameters: D\ = 17.4 nm, D cr = 1000 nm (BaTiC^); D\ = 8500 nm, D cr = 12 nm 
(PbTi0 3 ). 
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